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SUMMARY 

The problem of sound generation by vortices interacting with an arbitrary body in a low-speed flow has been 
investigated by the method of matched asymptotic expansions. For the purpose of this report, it is convenient to 
divide the problem into three parts. In the first part the mechanism of the vortex/body interaction, which is essen- 
tially the inner solution in the inner region, is examined. The trajectories for a system of vortices rotating about their 
centroid are found to undergo enormous changes after interaction; from this, some interesting properties emerged. In 
the second part, the problem is formulated, the outer solution is found, matching is implemented, and solutions for 
acoustic pressure are obtained. In the third part, Fourier integrals are evaluated and predicated results presented. An 
examination of these results reveals the following: (a) the background noise can be either augmented or attenuated 
by a body after interaction, (b) sound generated by vortex/body interaction obeys a scaling factor, (c) sound intensity 
can be reduced substantially by positioning the vortex system in the “favorable” side of the body instead of the “un- 
favorable” side, and (d) acoustic radiation from vortex/bluff-body interaction is less than that from vortex/airfoil 
interaction under most circumstances. 


1. INTRODUCTION 

Since the pioneer work of Lighthill, many theories and numerical methods for predicting aeroacoustics have 
been put forward. Most of these theories involve volume or surface integrals to account for source distributions in 
the turbulent near field. However, since a quantitative description of turbulent flows is difficult, perhaps we can, 
from a practical point of view, justify exploring some simpler methods. 

The aerodynamic aspect of sound production is known to be characterized by the formation of eddy motions 
(see Powell’s formulation of vortex sound in ref. 1 or 2). A simplified representation of these motions is the singular 
distribution of vorticity, such as the discrete vortices moving in an inviscid fluid. The justification for this can be 
found elsewhere (e.g., ref. 3); however, a brief description is provided here. 

In a viscous fluid and at a relatively low Reynolds number, the velocity and vorticity fields are usually continu- 
ous. As the Reynolds number increases, the vorticity field becomes more concentrated and is often embedded in an 
otherwise irrotational flow field. In taking the limit, these concentrated vorticity fields, which are governed by the 
Biot-Savort law, may then be modeled by discrete vortices with given circulations. The discovery of large-scale 
coherent structures in turbulent flows, especially in mixing layers, provides a further justification for this model. 
Considerable numerical work has been done in the application of various forms of discrete vortex methods to simu- 
late two-dimensional mixing layers (e.g., ref. 4). 

Because of a great increase in computing speed, aeroacoustics computations based on the unsteady Euler or 
Navier-Stokes equations can now be performed. However, obtaining accurate solutions depends on the resolution of 
certain issues, such as the size of the computational domain, the outer boundary conditions, and such. Thus far, no 
satisfactory solutions for these issues have been attained (ref. 5). 

Notwithstanding these problems, computational aeroacoustics may still be the most promising method in the 
long run. At present, the use of simpler methods without these problems may be justified. To this end, one chooses a 
potential flow model with discrete vortices to represent the flow in the inner region. And it is known that “all acous- 
tic motions in the vicinity of a singularity are solutions of Laplace’s equation” (p. 40, ref. 6). To solve acoustic prob- 
lems, the solution in the potential flow region has to be matched to the outer solution of the compressible flow 
(acoustic field) by means of the method of matched asymptotic expansion. The acoustic equation to be used herein is 
the acoustic equation for a moving medium instead of that for a stationary medium. 
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The problem now has been separated into two parts: an inner potential flow region and an outer acoustic field. 
This enables us to deal with complex geometries for the following two reasons, (a) There are a number of numerical 
methods available for two-dimensional potential flow in the presence of arbitrary bodies. The method used here is 
based on Martensen’s formulation and is that favored mostly by British researchers (e.g., see ref. 7 or 8). The advan- 
tage of this method for the present problem is that a body is replaced by a distribution of surface vorticity elements. 
Thus, the surface and the moving vortices can be treated similarly by the same Biot-Savart law. (b) The detailed 
flow characteristics of the inner region in an acoustically compact region will not, in general, affect aeroacoustics. 
Instead, only the gross properties of the potential flow influence the acoustic radiation. The evidence for this is that 
after a Fourier transform has been taken, the acoustic equation becomes a Bessel equation whose solution for the out- 
going wave is a Hankel function. Thus, the solution for sound radiation in the frequency domain is nearly determined. 
The matching process merely determines the coefficients. For this reason, even for a body of very complex geometry, 
once the inner solution is known, the pressure fluctuation in the far field in the frequency domain is established. 

The idea of using asymptotic matching to study aeroacoustics in a low-speed flow has been implemented in the 
past (e.g., refs. 9 to 15). In these studies if a body was present, it was usually of a simple shape, so the method of 
conformal mapping was applicable. Here we propose to broaden this method to include an arbitrary body and to 
adopt the acoustic equation in a moving medium. 

Although the outer solution is known in the frequency domain, an inverse Fourier transform is still necessary to 
bring it back to the physical domain. This has to be done numerically for a general case. There are at least two meth- 
ods available. One is the method of the discrete Fourier transform, and the other is the method of convolution inte- 
gral. The former is well known and its application to the present problem is relatively straightforward. The latter 
involves divergent integrals, which were derived by Conlisk and coworkers (ref. 1 3 and related papers), but which 
have not been evaluated numerically for far-field acoustic pressure fluctuations. Let us re-examine this method to 
see why divergent integrals appear, and to evaluate these integrals with the help of Hadamard’s (ref. 16) definition 
of the final part. 

This study examines sound generation by a vortex/body interaction. Initially a group of spinning vortices (see 
fig. 1) is introduced at an upstream station; as they move downstream and pass near the surface, an interaction takes 
place. This interaction causes the vortex trajectories to deviate from their original paths, alters the characteristics of 
the background noise, and gives rise to sound radiation. We studied each of these aspects in some detail and present 
the results herein. 

The results of this investigation are presented in the following sections. In sections 2 and 3, respectively, we 
derive the potential flow solutions in the inner region and discuss the related topics of vortex trajectories and their 
interaction with the body. In section 4, a brief description of the formulation is given, which leads to the Fourier 
transform of the inner solutions in section 5 and the outer solutions in the frequency domain in section 6. Section 7 
considers the process of matching and the solution for the acoustic pressure. Sections 8 to 10 cover two methods for 
evaluating the Fourier integrals and the test cases for the background noise. In section 1 1 results are presented and 
discussed. 


SYMBOLS 


c airfoil chord length 

ds m length of surface element at point s m 

K influence coefficient 

M total number of surface elements 

Af 0 free-stream Mach number 

N total number of free vortices 

p' fluctuating (acoustic) pressure 

R radial coordinate in outer region 
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z = x + iy 
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T 
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*F 
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CO 


radial coordinate in inner region 
position of surface element 
time in moving frame 
free-stream velocity 

velocity components in x and y directions 

surface velocity 

complex potential 

Cartesian coordinates 

complex number 

initial position of vortex centroid 

radius of spinning vortices 

strength of vortex 

surface vorticity per unit length at point s m 
radial coordinate of surface element s m 
angular coordinate of observation point 
angular coordinate in inner and outer regions 
angular coordinate of surface element m or n 
time in fixed frame 
velocity potential 
stream function 

angular velocity of spinning vortices 
angular frequency in Fourier transform 
frequency of spinning vortices 


Subscript 

i ith free vortex 

m,n surface element (point) m,n 

0 initial position or free-stream condition 
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2. POTENTIAL FLOW SOLUTIONS IN THE INNER REGION 


As stated previously, an incompressible potential flow is assumed to exist in the inner region, which includes a 
body as well as the moving rectilinear vortices with constant circulations concentrated at the singularities. The com- 
plex potential W for an isolated vortex located at zq a free space is 


W = <£ + /¥' = — log(z -Zo) (1) 

2 u 

where <p is the velocity potential; is the stream function; r is the circulation defined to be positive in the clock- 
wise direction; z = x + iy ; and z = x + vy . This is the building block for the potential flow to be discussed. 

Our present interest is to study the interaction of vortices with a body in a stream and, therefore, the presence of 
the body has to be taken into account. First, we consider Martensen’s method of potential flow calculation for a 
body of arbitrary shape without the moving vortices (for more details, see either of Martensen’s papers, or ref. 7 or 
8.) Since Martensen’s method is not based on a transformation method, no images of vortices are necessary. The 
body is defined by a string of surface vorticity elements separating the moving stream outside the body from a sta- 
tionary medium inside the body. The complex potential is simply a summation of M surface vorticity elements of 
undetermined strength, denoted as y m , plus the uniform free stream: 

M iy A* 

W = <P + i'F= yJ^-£Llog(z- Z|n ) + z (2 ) 

, 2 it 
m - 1 

Symbols here denote dimensionless quantities and are defined as follows: 

W = W/U 0 I, Z = z/i , Y m =rJu 0 , ds m =As m /l 

where Uq is the free-stream velocity; / is the reference length; y m is the surface vorticity per unit length located at 
the point z m on the surface of the body; and As m is the length of the element (see fig. 2). (In subsequent calcula- 
tions, the reference length for an airfoil is the chord and that for a rectangular block is the height.) The last term in 
equation (2) is the complex potential for the uniform free stream assumed to be parallel to the x-axis. It follows that 
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r m ds„ 


x-x„ 


m=l lK +(>’-3'm) 


( 3 ) 


The symbols u and v denote the velocity components in the x- and y-directions at a point x and y in the space in- 
duced by M surface vorticity elements. To determine the unknown y m in equation (2), we apply equation (3) at the 
surface elements and impose the Dirichlet condition that the internal tangential velocity component is zero: 
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( 4 ) 


Here x n , y n are the coordinates for R n is the radius of curvature at s n \ d n is the tangential angle at ^ (see fig. 2); 

and K(s n , s m ) is the influence coefficient, which represents the tangential velocity component induced at s n by other 
surface vorticity elements. The second term on the left side of the equation refers to the self-induced velocity of 
vortex element Y n ds n . 

Equation (4) is a system of M simultaneous equations for Af unknowns y m . The potential flow problem for the 
inner region is now reduced to the solution of M linear simultaneous equations. Any standard method can be used to 
solve it, for there is no ill-conditioning involved. We have used the Gaussian elimination method for every calcula- 
tion performed so far. 

The accuracy of this method has been assessed by Lewis (ref. 7), who compared its results with the known 
exact solutions, including those for bodies with sharp comers. We have also done our own comparisons, with such 
configurations as the NACA 0012 airfoil, and like Lewis, we concluded that we have a “high level of confidence” in 
this method. (For a cusped airfoil the accuracy may deteriorate.) Plotted in figure 3 is the velocity distribution at the 
forward comer of a sharp-edged rectangular block in a uniform flow. This is an indirect comparison since as the 
length increases, the potential flow at the sharp comer approaches that of a step for which the exact solution exists. 
The 2: 1 ratio rectangle can almost simulate a step. Also plotted in figure 3 is the surface velocity distribution in the 
forward face of the same rectangle but with round-edged comers. Both of these configurations will be used later for 
predicting sound generation, and we will see that even slight rounding of the comers can reduce sound radiation. 

Let us now consider the problem of vortex/body interaction. With the presence of free vortices in the field, 
equation (4) assumes a somewhat different form: 


2 n 


M 

X ymds m K{s„, 

m=\,m*n 


t) 4" Yn 


ds n 

4 7CR„ 


1 N 

= — COS Q n — — X, » Zi ) 


2 *f=l 


(5) 


The last term on the right side of this equation accounts for the presence of the N free vortices, whose influence co- 
efficient assumes the form 


w ^ (y n - y, ) cos fl n - (* n - x i ) sin e n 
/ \2 / \2 
(*n x i) 4 (y n )>i ) 


A two-dimensional vortex induces no velocity at its own center and moves with the local ffee-stream velocity 
plus the velocities induced by the other vortices. With the presence of these free vortices T), equation (3) becomes 


Uj-tVj 


M 

X 

m=] 


‘Ym ds m 


N 


in i 


2n 


Zj - z r 


+ 1+ X 

iJ7Uj 2n z J-n 


( 6 ) 


where z - is the coordinate for vortex j, and u jy Vj are the velocity components at z } - Here z i and y m are functions of 
time, but ds m , z m , and r are constants. 
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The evolution of vortex j is governed at each time step by the following two first-order differential equations: 



dyj 


dr 


(7) 


For a field with N vortices there are 2 N such equations whose solutions give the trajectories for jY vortices. At each 
time step, we have to update y m , by solving equation (5). A second-order method is used for computing the time 
development from t to x + Ax: 

Xj (t + At) = Xj (t) + J [3 Uj (t) - Uj( x - 4t)]4t 
yj{x + At) = yj (T) + J [3 vj (t) - v y - (r - At)]At 


The time step At is generally determined by comparing two solutions, one obtained with At and the other with 
At/ 2. If these two solutions agree reasonably well, the time step At is chosen; otherwise the process continues. In 
general, At= 0.0125 to 0.025 was found to be acceptable. The initial conditions for equation (7) are the coordinates 
of free vortices at locations sufficiently far upstream, where the interaction with the body is minimal. 

The accuracy of this method for a single free vortex was assessed by comparing two trajectories past a circular 
cylinder in an incompressible flow; one was obtained by the present method and the other was the exact solution 
obtained by the image method. In a test run with At = 0.1 and T= -1.0, four-figure accuracy was attained in the 
vortex position, except very near the surface of the body, where the accuracy was reduced to three figures. 


3. VORTEX TRAJECTORIES AND INTERACTIONS 

With the necessary equations for a potential flow in place, we can now proceed with the solutions. First, a few 
remarks are in order. 

Three conservation laws apply to a system of discrete vortices in a two-dimensional unbounded fluid at rest at 
infinity; one of these states that the centroid of the system is stationary. If there are N equal, like-signed vortices 
spaced equally around the circumference of a circle, they rotate about the centroid with a constant angular velocity. 
In the present study N is equal to 4 or less. Emphasis will, however, be placed on a pair of spinning vortices, since 
this is a classical example for acoustic radiation and there is an analytic solution for it. When such a system is placed 
in a moving stream, the centroid moves with the local convective velocity. 

The central topics of the present investigation are the vortex trajectory patterns and the acoustic radiation from 
the interaction. Problems of a similar nature have been previously studied. For example, Hardin and Mason (ref. 17) 
examined the blade/vortex interaction for a single vortex and the associated sound generation. Their results will be 
used later for comparison with the present solutions. More recently, Panaras (ref. 18) investigated the vortex trajec- 
tories for many vortices and their interactions with several Joukowski airfoils. However, the relatively large number 
of vortices used in the investigation may have obscured the process, since the vortices developed irregular patterns 
prior to the interaction. As a result, the effect of the interaction on the vortex trajectories cannot be easily seen. 

In order to see the effect of an interaction, it is preferable to limit the vortices to a small number and to arrange 
them in such a way that the motion prior to interaction is regular; thus, the change in their trajectories can be dis- 
cerned immediately. To this end, the present system consists of no more than four vortices, and it undergoes a circu- 
lar motion in a stationary medium. With the convective velocity included, the trajectory of each individual vortex 
forms a wave pattern (figs. 4 to 9). 

Hardin and Mason pointed out that vortex/blade interactions are highly nonlinear, which is also true here. So a 
slight change in the initial conditions may result in distinctly different trajectories after interaction. For this reason a 
number of plots are included here to show the different possibilities that an interaction can bring. In figure 4 three 
trajectories by a pair of spinning vortices of equal and positive strength (clockwise rotation) are shown. Initially this 
pair was stacked in an upstream station (fig. 1), whose vertical distances measured from the airfoil chord line are the 
same, but the horizontal distances from the leading edge of the airfoil are different. With jcj = x 2 = -24.4, the fore 
and aft trajectories are similar. However, as the initial position of the pair moves forward to -24.6 (20 percent of the 
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chord), the trajectories after interaction (the aft trajectories) are no longer the same as the upstream trajectories (the 
fore trajectories). The pair spins faster behind the airfoil. If the pair moves further upstream to -25.0, a pronounced 
change in the aft trajectory pattern takes place, in which case the vortices depart from each other and, consequently, 
the rotating speed becomes much smaller. This suggests that the background noise of a spinning vortex pair can be 
attenuated by interaction with a body, which we will examine later. To understand this scenario, we cite a well- 
known equation given in reference 19: 


i2 = — J =T (JV-l)r (8) 

Am 

where Q is the angular velocity and N represents N equal, like- signed vortices r arranged uniformly about a circle 
of radius a . According to this equation, the spin rate becomes smaller as the two vortices depart from each other, 
and the associated sound decreases. 

A closer examination of the trajectories reveals that here the aft trajectory patterns are periodic in the sense that 
their evolution from a high to a low rotating rate repeats itself as the initial position of the vortex pair is repositioned 
continuously forward or backward. The reason for this is that a vortex moves along a wavy path and its time of 
arrival in the neighborhood of the airfoil is periodic. 

Thus far, only the positive vortices placed above the chord line have been studied. If they are replaced by two 
negative (counterclockwise) vortices with the other conditions unchanged, the fore and aft trajectories of all three 
previously mentioned cases become similar. For example, figure 5 is a counterpart of figure 4(c), with identical 
parameters except for the sense of the vortices. In reference 14, Sen quotes Howe’s result that a vortex produces 
sound mainly when it cuts across base streamlines near an edge. The velocity field induced by a negative vortex is in 
the direction of the local (base) convective velocity and, therefore, produces less disturbance. By contrast, the veloc- 
ity field induced by a positive vortex goes against the local convective velocity near the surface, so the vortices are 
more likely to be repelled to a different orbit. On the basis of this intuitive notion, we would expect a pair of positive 
vortices placed below the chord line to have the same effect as a pair of negative vortices above the chord line. This 
expectation turned out to be correct and has been verified by calculations. 

Three different trajectory patterns are shown in figure 6, one of which involves four vortices. They are included 
here because these patterns will be used later to calculate acoustic pressures. Note that the bluff body (a 2: 1 ratio 
block) shown in figure 6(b) and the following figures is merely a distorted schematic of the real geometry. Since 
there is a large disparity between the vertical and horizontal scale, a 2: 1 block would be unrecognizable if plotted to 
scale. 

Apparently this intuitive notion for an airfoil is also applicable to the vortex/bluff-body interaction. For 
example, see figure 7, where two sets of trajectories for a sharp-edged bluff body are plotted. Again, the initial con- 
ditions for parts (a) and (b) of figure 7 are the same; only the sense of the vortices is different. 

Previously, we dealt exclusively with equal and like-signed vortices. If we ease this constraint a little to permit 
vortices of unequal strengths, the centroid will be shifted to one side; other aspects of the trajectory pattern will, 
however, remain essentially the same. For this reason, no results and graphics have been included. 

Let us now shift our attention to the interactions involving three like-signed vortices of equal strength rotating at 
a constant speed (fig. 1). In this case the forward trajectories are essentially three intertwined identical curves, but 
the aft trajectories are a variety of vastly different patterns. Consequently, only a representative case is shown in 
figure 8, which consists of two sets of curves, one for a sharp-edged body and the other for an equivalent round- 
edged body. The purpose is to demonstrate the influence of a comer. Apparently, rounding off the comers slightly 
can make a difference in the aft pattern as well as in the acoustic radiation. 

Usually, after a vortex/body interaction, the trajectories become less organized, as can be seen in the aforemen- 
tioned figures. The reverse may also be true in the sense that after an interaction vortices may spin in a more orderly 
fashion than before. This possibility is depicted in figure 9, which shows all three aft vortices rotating at nearly the 
same speed. 


4. FORMULATION 

The problem to be considered consists of two regions: an inner and an outer region. The governing equations 
for the former have been given. In this section the differential equations for the latter are to be derived. Since the 
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acoustic pressure radiated from a source will be distorted after traveling through a convective stream, it is preferable 
to find solutions based on the acoustic equation in a moving medium instead of in the stationary medium used in 
earlier investigations. This equation for a fixed frame of reference is given by 


<?V 2 P 2 P' 

dx 2 ° Dx 2 


0, where 


_D_ = d_ ± 

Dz dx dx 


( 9 ) 


The symbols p\ T, x y y, Af 0 are, respectively, the dimensionless quantities for the perturbed pressure, time, spa- 
tial coordinates, and Mach number. They are defined as follows: 


P’ = ^=T 

PoVo 


U 0 t x y %a U 0 

t = -=K X ~T’ y = T' M o=- 

l l l c 0 


where p 0 is the density, l is the reference length, / is time, and c 0 is the speed of sound. To make the analysis more 
tractable, this equation is transformed from a fixed frame to a moving coordinate system (ref. 20) by defining 
x = x - T, y' = y, and t = x. Thus, 


d 2 i 


d 2 i 


dx' 2 dy' 


.2 


-Mi 


2 d 2 p' 
dt 2 


= 0 


( 10 ) 


This equation is of the same form as the classical acoustic equation. Therefore, all available methods for solving the 
classical equation can now be put to use. Once the solution is known, it can be transformed back to the fixed frame 
of reference. 

Introducing two outer variables defined as R = M 0 r' and 9 ' = tan -1 y'/x\ where r' = ( x ' 2 + y' 2 ) 1/2 , we obtain 


d 2 p' 1 d p' 1 d 2 p' d 2 p' 

jR 2+ l~dR + ^Jff Y ~d7 I 


( 11 ) 


This is the final form of the governing equation for the outer region, subject to the radiation condition at infinity and 
asymptotic matching near the body. 

On the basis of the premise that steady flows generate no sound, the inner potential flow equations and the sub- 
sequent numerical expressions can be written in such a way that the unsteady part is separated from the steady part. 
For the purpose of matching, we have to express the unsteady velocity potential in terms of p' given by the linear- 
ized Bernoulli equation: 


P 


= P~Po =~ 


D&' 

Dx 


d&' 

dt 


( 12 ) 


The operator D/Dx represents the time rate of change of (P'as seen by an observer moving with the free-stream 
velocity. Note that the notation x denotes time in a fixed frame. If a moving frame is used, Twill become t and 
D&'/Dx becomes d&'/dt. 


5. FOURIER TRANSFORM OF INNER SOLUTION 

With the governing equations for the inner and outer regions both attained, we can proceed to solve low- speed 
aeroacoustics problems. As briefly mentioned in section 1, the procedure is as follows. First the unsteady part of the 
inner solution in the limit as r — » °o is determined, from which, by means of equation (12), the perturbation pressure 
is found. Next, the Fourier transform of p ' is obtained for matching. The partial differential equation (eq. (11)), is 
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then transformed via the Fourier transform to an equivalent ordinary differential (Bessel) equation, whose solution 
for the outgoing wave at infinity in the frequency domain is the Hankel function. The limiting form of this function 
as R — » 0 is matched to the inner solution to secure the unknown coefficients. Then in the last step, an inverse trans- 
form is performed to return the solution to the physical space. 

Following these steps, let us proceed to determine the inner solution by first considering a pair of spinning vor- 
tices moving past a body. The complex potential for this system is similar to that of equation (2): 


W = <P + W=Y j !^£ SL \og(z-Z m ) + Z + lf-log[z-(T + Z cl +Zo+Zij\ + !f-\og[z-{t + Z cl -Zo + Z2)] (13) 
, 2k 2 K IK 

m - 1 

Most symbols in this equation have been previously defined. The terms in the last two sets of brackets represent the 
positions of the spinning pair. The symbol z ct refers to the initial position of the vortex centroid and is a constant. 

The quantity T + z ct is the position of the centroid at time r converted by the free stream in the absence of a body. 

The symbols ±z 0 are the local coordinates of vortices relative to the centroid (fig. 1). If the body is absent, this pair 
rotate freely, are always aligned, and need only the quantities r + z ct ± Zq to denote their instantaneous positions. 
However, in the presence of a body, these positions will be perturbed because of the interaction with the body; thus 
we introduce Zj and z 2 

After excluding thie steady part in equation (13), the relevant portion is the unsteady part W given by 


W' 


M 

■S 

m~ 1 


l lm^m |Qg | 
2 K 


z\\-^ 

Z 


iT_ 

In 


log 1 - 


Z-T-Z ct 


+ log 1 + 


ZQ-Z 2 

Z-T- 


z ct ) j 


Expanding this equation for I z^/z I « 1, 1 (Zq + z{)!{z — T— z ct ) I « 1*1 (z 0 -z 2 )/(z - t- z ct ) I « 1, and defining z - r 
exp (iO) and z m = S m exp we have 


W'- 


M 

X 

m - 1 


2k 


log r + id - — ^ e‘^ m ^ 
r 


2k 


*2 


zl 


(Z-T-Z c< ) (Z-T-Z cr ) ( z - T - 2cr ) 


+ ... (14) 


The first two terms in the first pair of brackets are zero, because V M Ym^ s m ~ ' s due to the fact that 

^ M y m ds m = ^ (d 0 / d s)ds = 0, provided that there is no circulation. 1 ( y m can also denote the surface tangential 
velocity component; see ref. 7.) Note that only one second-order term in this equation is retained; others have been 
dropped. This second-order term represents the unsteady motion of the spin and produces the background noise. 
Other terms in this equation are first-order terms that are essentially zero until the vortices experience an interaction 
with the body, and they will return to zero after the interaction is over. 

The next step is to put equation (14) in a more convenient form and then to find// by means of equation (12). 
Introducing the polar coordinates of the following quantities, 

z - r - z ct - f 3 e l ° and Zq = P e l<0Qt 

where co 0 = I747t P 2 , we transform equation (14) to the form 

= 0 , ■ ■ y Ym ds m -e) »T x x + iy, + x 2 + iyj e -iff _ir_ff_ e li(<o 0 ,-6-) (1 , 

" 2k r 2 k r' 2 k r ’ 2 

m~ 1 


l If there is a strong vortex system in the proximity of a symmetric airfoil, the airfoil may experience a local angle of attack even if the free stream 
is parallel to the chord line. This local effect, however, is neglected here. 
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In this equation the symbols r and Ain the first term are the polar coordinates measured from the fixed frame, since 
the body is stationary. However, the acoustic equation in the outer region is given in terms of / and ff. Thus, it is 
necessary to express r and Ain terms of /, ff, and t before any matching can take place. Although this is possible, 
we can demonstrate mathematically to first-order accuracy that r and A can be replaced by / and ff during the pro- 
cess of matching. The physical explanation for this is that dyjdt ± 0 only when vortices are in the vicinity of the 
body, where r and flare approximately equal to r and ff. Note also that by defining = /? exp (icotf) we implicitly 

assume a semblance between the fore and aft trajectory patterns; otherwise z 0 has to be defined differently. 

Taking the real part and applying equation (12), we obtain 


M 


Pi = 


m - 1 


jYm ds m §n 
d T 2 n r' 


sin(0'-0 m ) + 


r 

2k r' 


Dx x 

~Dr 


Dx 2 
Dt 


1 sin O' - cos 6' - — -L- cos 2(<u 0 f _ 8')- 

) \ Dt Dt ) J K r' 2 K ’ 


(16) 


The first term represents the fluctuation of the surface pressure, the second term is due to the interaction, and the last 
term represents the unsteady motion of the spinning pair. Notice that y m is differentiated with respect to r since the 
body is fixed. Applying the Fourier transform pair 


/(/?,<»)= jf(R,t)e~ i0X dt and f(R,t) = ^- jf(^ 0 ))e ,OJ, dco 


to equation (16) gives the Fourier transform of p-’: 


> _ Mo V ds m d m 
Pi 2 k ^ R 

m= 1 


sin 






K 


Qy\ + Dyj 1 „-ica 

Dr Dr 


*dt 


dt 

q2 


2 71 


sin A' r 

Dx x 

Dx 2 
+ — - 


R J l 

— CO 

,Dt 

D T, 

1 R 


- | cos 2 $'^S(o) _ 2co 0 ) + 8(g) + 2g) 0 )] + i sin 2A'[5(g) + 2cu 0 ) - 8(co - 2g) q )]J. 


(17) 


where R = Af^'and 5 ( ) refers to the 5-function. The second-order term is proportional to MR 2 , whereas the first- 
order term is proportional to MR. 


6. OUTER SOLUTIONS IN THE FREQUENCY DOMAIN 
In view of the inner solution (eq. (17)), the outer solution may be assumed to have the following form: 

M 

p' 0 (R,6',t) = ^ (/?, t) sin(fl' - fl m ) + /> 1 (/?,f)sinfl' + p 2 (Rj) cos 8' + pi(R,t) sin 26' + p 4 (R,t)cos26' (18) 

m - 1 


Substituting this expression into equation (11) gives the following equations: 


d\_ + }_dq L 

dR 2 R dR 


1 _l3L 

R 2 ^ dt 1 


= 0, 


d 2 


?2 + J_ 


iq 2 


dR 2 RdR R 


2*2 


^ = 0 

dt 2 


where q j refers to p m , p x , or p 2 , and q 2 refers to /? 3 or p 4 . These two equations can be transformed via the Fourier 
transform to two equivalent Bessel equations: 
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d 2 q x 1 dq { ( 2 1 

dR 2 R dR l, R‘ 


(O <?i = 0 and 4- + — 

R 2 ) dR 2 R 




whose solutions for the outgoing waves at infinity are the Hankel functions of the second kind. We, therefore, obtain 

Po = ^A m Hl 2 \caR)sin(6' -6 m ) + AqH^^coR) sin O' + Bo//, (2) (ftV?)cos0' + C 0 H (2 \caR)sm 20' 

m=1 +D 0 H^ ) (aR)cos26' (19) 

where A m , A 0 , ... D 0 are the unknown coefficients to be determined by matching. 

7. ASYMPTOTIC MATCHING 

As mentioned briefly in Section 5, matching requires that the outer solution be put in the limiting form as R — » 0 
(ref. 21). With the aid of the formulas given in reference 22, we have 

M ry • 9 / 4/ 

Pq ~ V — — A m sin(0' - 6 m ) + ( Aq sin 6' + Bq cos 6') + — j{Co sin 20* + A) cos 2 O') 

F0 ■“ kcoR m v ’ ncaR v ^ ’ k(d 2 R 2 v 

m=l 

Equating this equation to equation (17) gives the unknown coefficients in equation (19): 


A/n = ~ 


iM o© * f dy m - ia ». 

— d*nfi m C 


irM^co 7 ( Dx l | Dx 2 \-ot 
^ 4 J \ Dr Dx) 


Bn = ~ 


im Q 0) r f Dy 1 


\--ia 
Dx f 


C 0 = — 0)0 to 2 [S(o) + 2© 0 ) - 5(o> - 20)0 )] 


+ i^ + j_iV _iV=o (11) 

R 2 ^ff 2 

With the coefficients defined, we apply the Fourier inverse transform to equation (19) to give the acoustic pres- 
sure in the far field: 

p, _ _<Mo ^ds m 8 m s j n (0' - e m ) J adi^ (aR)e ia * dG> J ^~e~‘°^d^ - sinfl' J otf/p^OlRJe'^dO) 

m-1 — oo — OO ^ 

x 7f^l+^2l + i™0 cosr ]o<>(otfy“</o> + 

+ + ®o') _ Y 2 (2o} 0 R)co&2(6’ + oty)] (21) 
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The last term, written in terms of Bessel functions of the first and the second kind, 7 2 (26>^) and Y 1 {2co ( fi), 
refers to the background noise that will be heard by the observer constantly. All other terms represent the interaction 
sound, which will be heard by the observer in the interval when the vortices are near the body. 

The remaining task is to evaluate the Fourier integrals in equation (21). Two procedures are generally available. 
One is to use the discrete Fourier transform, and the other is to rewrite the integrals as convolution integrals. These 
convolution integrals are akin to the divergent integrals of the Hadamard type and have to be evaluated accordingly. 
These two procedures will be discussed separately. 


8. METHOD OF DISCRETE FOURIER TRANSFORM 

Although the term coH^ 2 \q}R) in the inverse integrals in equation (21) approaches co l/2 as co — » «>, the pulse 
functions such as dyjd^ are basically of the triangular wave forms whose Fourier transforms decay like 1/m 2 at 
infinity. Thus, the function as a whole is absolutely integrable, and the inverse integral exists. On this basis, the 
existing methods for the discrete Fourier transform, such as those stated in reference 23, can be used. In order to 
provide some assurance that this method was being applied properly, we repeated Crighton’s solution numerically 
(ref. 10) and found a close agreement with his results. 


9. METHOD OF CONVOLUTION INTEGRAL 

The procedure of writing the Fourier integrals in equation (21) in terms of convolution integrals was developed 
by Conlisk and coworkers (see ref. 13 and related papers). To our knowledge, no actual evaluation of these integrals 
has been made. Their concern has been with the near-field sound, which does not require evaluation. A close look at 
their procedure suggests some questionable steps may have been taken. Such steps could be the reason that the 
resulting integrals become divergent and are of the Hadamard type. However, the evidence that we have, on the 
basis of actual evaluation, seems to indicate that these steps did not totally invalidate the outcome, except the emer- 
gence of divergent integrals. With the help of Hadamard’ s definition of the final part (ref. 16 or 24), the physically 
meaningful part of the integral is extracted and the meaningless part is discarded. Perhaps for this reason, the acous- 
tic pressure in the far field based on convolution integrals agrees reasonably well with some existing solutions under 
certain circumstances. However, discrepancies begin to develop if the interaction becomes strong and vortices are 
very close to the surface. This anomaly may be common to the methods using discrete vortex filaments with singular 
centers, as was also mentioned by Hardin and Mason. 

In order to make the explanation easier, let us repeat some of the steps in reference 13. Equation (21) contains 
many terms, but they are all similar. It is, therefore, sufficient to pick only one term for illustration: 




^ajH {2 \coR)e l<a d(0 J ~^ € 


7 ( 22 ) 

8tt J 1 

— oo — oo 


The purpose of this step is to change the order of integration, which is permissible if the functions are square - 
integrable (see p. 27 in ref. 25). However, as mentioned previously, coH^ 2 \coR) — » o^ 12 as Gt) — > The square- 

integrability criterion do)<°° i s therefore not satisfied. 

— oo 

Assuming this step is valid, we still have to write H^ 2 \o)R) in terms of the modified Bessel function K } ( ) be- 
fore we can consult a table. Equation (22) can now be reduced to 


\wH[ 2 \<oR)e i(a d(o] = ^d^\(0K,{i(0R)e ,coU ^da 


m 0 R ‘ j R Dx x 


dt 


4n 




(23) 
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Although we formally obtain a convolution integral, an improper use of the table has been made. For example, refer- 
ence 26 states explicitly that the parameter in AT,( ) must be real and positive. This constraint is violated in the above 
equation. We nevertheless consider that this method has its merit as will be demonstrated later. 

The convolution integral in equation (23) is similar to the acoustic potential function given in reference 13, 
which has to be differentiated in order to obtain the acoustic pressure. However, special care is needed when a diver- 
gent integral is involved. For this reason p' is chosen as the dependent variable instead of the acoustic potential. 

Divergent integrals similar to that in equation (23) have also been encountered in the linearized theory of super- 
sonic flow (slender-body theory), which is known to be equivalent to the acoustic theory. Furthermore, the present 
convolution integral is of the Hadamard type. Therefore, existing methods (ref. 16 or 24, for example) can be used to 
define the final part. A brief description is given here to better illustrate a point that may affect the accuracy of some 
numerical results. 

Using/(^) to represent Dxj/D£, we may write equation (23) as 


f " r 

Jr . o 9i 3/2 


t-R 


/($-/(*-*) 
3/2 


t-R 




dS + f(t-R) J 


dS 


r o ,i3/2 
[( t - S ) 2 -/? 2 


The first integral on the right side of the equation can be shown to be an ordinary improper integral and may be 
evaluated accordingly. The second integral is a divergent integral and has to be determined by using Hadamard’ s 
definition of the final part. We now write the first integral to be a sum of two parts with o denoting a small quantity 


t-R 

j 


m-fit-K) 

3/2 a( => 


t-R-c 


[«-tf-R 2 ] 


m-AL-R). „ 

i3/2 


t-R 


AS) -f (t-R) 

3/2 




(24) 


The first part is a regular integral and can be evaluated by a conventional method; the second part, after expanding 
f(S) about the point t-R and retaining the first-order term, becomes 


J i ?l3/2 

t-R-c (t~&) —R 


t-R 

d% ~ ~f'(t - R) J 

t-R-a 


it-R-S) 

(t-£-R) V2 (t-£ + R) y2 


d4 


which has a closed-form solution. Our experience seems to indicate that crin equation (24) has to be chosen in such 
a way that it is comparable to the size of Aij. If it is too small, large local spikes may develop. Even with some pre- 
caution, this anomaly is still there and can be seen in figure 10. 


10. BACKGROUND NOISE 

It is perhaps worthwhile at this stage to present some evidence that the method of convolution integrals is 
capable of yielding reasonable results, at least under some conditions. To this end a pair of spinning vortices con- 
vened in a free space is chosen as the first test case. An analytic solution to this problem is known (ref. 9) and is, in 
fact, the last term in equation (21). The process leading to an analytic solution involves the use of the 5- function. If 
we forgo this practice, we can still proceed, but in this case convolution integrals will emerge. Taking similar steps 
as in equations (22) to (24), we have 


:~M 2 R 2 co 0 li 2 
4 71 


t-R 

sin 28' J 


sin 2 0)qS 


t-R 


[(r-£) 2 -tf 2 ] 5/2 


d£ + cos20' J 


cos 2 (O 0 S 
[(r-£) 2 -/? 2 ] 


<%\ 


(25) 
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This expression is similar to equation (23), except that the denominator in the integrand is raised to 5/2 power 
instead of 3/2 power, which occurs because the acoustic pressure for a spinning pair is a second-order effect. 

Figure 10 compares two sets of p . Although the frequency and the amplitude given by equation (25) are both 
somewhat larger than the analytic quantities, the agreement is not really bad, especially in view of the fact that the 
solution given in equation (25) involves two divergent integrals. In the present formulation, which uses the convec- 
tive acoustic equation, the Doppler effect on acoustic radiation should appear by itself. Although this property can- 
not be seen easily in figure 10, it will be discernible at a higher Mach number — for example, at Af 0 > 0.3. 

The acoustic pressure depicted in figure 10 may be assumed to represent the background noise in the problem of 
a vortex/body interaction, if the fore and aft trajectories are similar. If they are different, the characteristics of the 
background noise after interaction will be altered, that is, have a different frequency and amplitude. In this case we 
have to choose two sets of average values to approximate the fore pattern from -<» to Xj and the aft pattern from 
x 2 to <». (The interval between x x and x 2 is the interaction zone, which is neglected as far as the background noise is 
concerned.) For illustration, the step modulated input /r) = U(t) cos co^t is used to represent a cosine wave in the 
interval from 0 to <*>, where U(t) is the unit step function. The Fourier transform of fit) is 

k<o) = j [S(a> - a> 0 ) + + “>0 )] + 

(See p. 41 of ref. 25.) This expression is seen to be equivalent to D 0 in equation (20). However, since it has a singu- 
lar term, a special treatment is required in the inverse transform. For this reason, the method of convolution integrals 
is preferred. 

Again following similar steps as in equations (22) to (24), we obtain 


P =- 


3fMp/? 2 

4n 


sin 20' 


2 f*t/(-g-.Xi)sin 2 


<»lPl J 


i2 ' r R U(^-x 2 )sin2co^ 


[(t-tf-R 2 } 5 ' 2 


d £ + (OjPi j 


[(r-S) 2 -* 2 f 




, 3I~MnR 2 

y + — cos 20 

4k 




t~Rr 


[it-tf-R 2 ] 512 




(26) 


where co x and p x denote the average angular velocity and the radius of rotation in the forward portion of the trajec- 
tory, and CO} and /? 3 denote the same for the aft portion. The symbol -x x is the endpoint of the forward portion, 
whereas is the starting point of the aft portion (see fig. 4(c)). The symbol U ( ) refers to the unit step function 
shifted by the amount -x x or x 2 to account for the semi-infinite length. If the segment becomes infinitely long, the 
unit step function becomes unity throughout and equation (26) reduces to equation (25). 

To demonstrate the applicability of equation (26), we extricated the required information from the calculated 
trajectories in figure 4(c) to evaluate p\ which is plotted in figure 1 1 . Notice that in this case the background noise 
after an interaction has been almost totally attenuated. Because of the linearity property, this noise signature can be 
superimposed on the acoustic radiation from an interaction to present a complete acoustic field. This is shown in the 
following section. However, because of the scale difference, the background noise is not discernible. For this reason, 
the body in figure 1 1 is removed and the vertical scale is enlarged. 


11. RESULTS AND DISCUSSION 

In the foregoing sections we discussed two methods to evaluate the Fourier integrals. The method of convolu- 
tion integrals is rather complex and lacks mathematical rigor, yet it is an interesting approach and may provide sim- 
pler solutions on occasion, such as in the previous example. Furthermore, two more cases will be given to 
demonstrate its applicability to sound generation, provided that the interaction is moderate. For a strong interaction 
overpredictions are liable to occur. Thus, the discrete Fourier transform (DFT) is our primary method. 
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As mentioned earlier, Hardin and Mason investigated the problem of sound generation by a vortex and pre- 
sented their calculated results for a Joukowski airfoil. Using the initial and boundary conditions indicated in refer- 
ence 17, we simplified equation (21) for a single vortex and solved it by means of the DFT to secure a set of 
corresponding acoustic pressures. Figure 12 shows this comparison. The characteristic length in reference 17 is the 
span length. However, no such length is available in the present investigation. Consequently, only the relative mag- 
nitudes can be compared. The overall shape appears to be similar, though details are different. Although Hardin and 
Mason analyzed both lift and nonlifting airfoils, only the nonlifting case is compared, since the angle of attack is 
assumed to be zero in this study. 

Observing the trajectories in figure 4, we could assume that the degree of interaction becomes progressively 
stronger from parts (a) to (c) of figure 4. We intend to show that as far as sound generation is concerned, this as- 
sumption is correct. For this purpose we used the DFT to solve equation (21). There is, however, a disparity in the 
fore and aft trajectories in parts (b) and (c) of figure 4; the last term in equation (21) for these two cases has to be 
replaced by equation (26). The results are plotted in figure 13, which shows that the background noise is relatively 
low, especially in figure 13(c). To discern the background noise in figure 13(c), refer to figure 11, which has been 
magnified. These three examples illustrate that the background noise is strongly affected by the dynamics of interac- 
tion, which can augment or attenuate its intensity. In the case of figure 13(c), the interaction almost silenced it. 

Let us now digress to test the accuracy of the method of convolution integrals by repeating the calculation for 
the case in figure 13(a). The results are shown in figure 14(a). The characteristics of the curves in these two figures 
are approximately the same, although the intensities are dissimilar to various degrees. We take a more lenient view 
of this method and consider this kind of accordance to be satisfactory. 

Returning to the method of DFT, let us consider the case in figure 5. The initial and boundary conditions for this 
case are identical to those in figure 13(c), except that r i = r 2 = - 0.15. Simply changing the sign of the vortices 
greatly diminishes the pressure fluctuation in figure 14(b), and the first peak is now positive instead of negative. As 
mentioned earlier, this outcome is borne out by Howe’s result. In other words, sound generation by vortex/body 
interaction depends on the sense of rotation. The effect is antisymmetric with respect to the chord line. For example, 
if the vortex in figure 5 is shifted to a corresponding position below the chord line, the trajectory will be the same as 
in figure 4(c), and the sound radiation will be identical to that in figure 13(c). 

A possible application of this property is that if a design change is feasible, we can position the vortices on the 
favorable side of the body to substantially reduce the noise level. Even for a single vortex, the direction of rotation 
can make some difference. For example, the predictions by Hardin and Mason in figure 5 of reference 17 show the 
different levels of sound radiation for Y 0 f c = -0.2 and +0.2, even at zero angle of attack. 

A pair of fast-spinning vortices would be expected to produce a higher level of sound radiation than would a 
pair of slow-spinning vortices. In order to substantiate this assumption, the initial positions of the vortices in figure 
4(b) was modified slightly by reducing the radius of rotation, thereby making them spin faster, as depicted in figure 
6(a). The corresponding acoustic radiation calculated by the DFT is plotted in figure 15(a); it is considerably higher 
than that in figure 13(b). 

Until now, the point of observation was assumed to be situated at 0= 75° from the positive *-axis (the chord 
line behind the airfoil). It would be interesting to know the sound radiation received by an observer at the same dis- 
tance from the center but who is located ahead of the airfoil and slightly above the chord line (0= 170°). In this 
case, as displayed in figure 15(b), the intensity is reduced in comparison to figure 13(c), the signals are bunched 
together, and there is a time delay in receiving the signal. Such characteristics are what we would expect. 

Let us now direct our attention to bluff bodies. Contrary to our intuition, it seems that sound radiation from 
vortex/bluff-body interactions is almost always lower than that from the vortex/airfoil interaction. The explanation 
lies perhaps in the flow pattern in the nose region of a bluff body, where the fluid, after impacting a flat surface, 
branches off into two streams and skirts the comers to be away from the body. As a result, the vortex never gets a 
chance to be very close to the surface, and there is less interaction. This seems to be true even for the almost head-on 
collision shown in figure 7(a). By contrast, a pair of spinning vortices can come very near the airfoil and practically 
“scrub” the surface, thereby giving rise to a strong interaction. However, it is uncertain whether the discrete vortex 
method is still valid, as vortices move into the boundary layer and the viscous effect becomes important. Further 
investigations may be needed. With an airfoil, the trailing edge may also intensify the interaction. 

In order to clarify some of these intuitive explanations, in figure 16(a) we plotted acoustic pressure fluctuations 
for the trajectory depicted in figure 7(a). As can be seen, this pattern changed enormously after interaction. In the 
case for the airfoil, such a pattern would usually indicate high intensity radiation. Yet the acoustic pressure in 
figure 16(a) is only moderate. (Note here that the background noise is not plotted in this figure and will not be plot- 
ted hereafter since its absence will not hinder our future discussion.) 
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As a further test of the method of convolution integral, another example is provided. This method was used to 
recalculate the acoustic pressure in figure 16(a). The results (see fig. 16(b)) show that the details are again at vari- 
ance with figure 16(a), but the intensity is comparable. 

As in figure 13, we can increase the intensity of the acoustic pressure by manipulating the vortex pair's horizon- 
tal distance from the leading edge. With a proper shift, the acoustic pressure can be increased considerably (see 
fig. 17(a)), but it is still significantly less than that in figure 13(c). To lower the intensity further, simply round off 
the comers slightly as in figure 3(b) or figure 8(b), and the pressure fluctuation (fig. 17(b)) lessens substantially. By 
shifting the vortex pair somewhat upward and away from the centerline, the serrate peaks in figure 17 give way to 
the smooth curve in figure 18. This transition from a serrate configuration to a smooth curve also occurs in a vortex/ 
airfoil interaction. 

Finally, let us investigate the cases of three and four spinning vortices. In order to explain the scaling effect of 
interaction, it is more convenient to pick the system of four vortices (fig. 6(c)). As in the previous cases, these are 
equal and like-signed vortices that are spaced uniformly around a circle (fig. 1 ) and spin at a constant rate in a sta- 
tionary medium (eq. (8)). Now we can remove two trailing vortices to form a pair of stacked vortices, as in previous 
dispositions, thereby making the rotation radius smaller. The acoustic radiation from the interaction of this pair with 
the airfoil is plotted as a solid line in figure 19(b). Next, the two previously removed vortices are restored to form a 
quadruple. This system, however, spins at a rate 1 .5 times faster than the rate of the pair, so the acoustic pressure 
becomes substantially higher as shown in figure 19(b). (The motion of four vortices in the present configuration is 
known to be stable (ref. 19), and there is no computational evidence of instability after interaction.) 

It is proposed that the increase in acoustic pressure with the quadruple (fig. 19(b)) is almost entirely due to the 
increase in rotational speed and nearly independent of the number of vortices involved. To illustrate this seeding 
factor, let us reduce the vortex strength for the quadruple from 0.15 to 0.1, so that the rate of spinning matches that 
of the pair, and recompute the acoustic pressure. This is represented by the dashed line in figure 19(a), which differs 
only moderately from the solid line. This scaling effect should also be valid for a system of three vortices and for 
bodies other than an airfoil; both of these cases are demonstrated next. 

Knowing the trajectory pattern in figure 6(b) and following a procedure similar to that just discussed, we plotted 
the solid lines in figure 20(a) to represent the acoustic pressure after the third trailing vortex is removed (fig. 1). 
Restoring this vortex to the system produces radiation that is much more intense (see fig. 20(b)). If we then reduce 
the vortex strength in a similar manner to make the two rotating speeds identical, we obtain the dashed lines shown 
in figure 20(a). These two sets of lines agree fairly well, so the scaling factor is again validated. 

It is pertinent to make another remark here. The scaling factor is expected to hold true for bluff bodies under 
most circumstances; however, it may not be valid for strong vortex/airfoil interactions because one or more vortices 
may scrub the surface. If the number of vortices is reduced, the likelihood of a close encounter is different. 

Although some indirect comparisons have been made herein, there is a genuine need for detailed comparisons 
with experimental data. The challenge to developing an experiment that simulates a two-dimensional flow may be 
alleviated if three-dimensional problems can be solved. This is feasible, especially if the body is relatively simple 
and the self-induced velocity on the curved vortex filament is calculated approximately. See, for example, 
reference 27. 


12. CONCLUDING REMARKS 

A procedure to solve aeroacoustics problems due to vortex/body interactions in low-speed flow has been pre- 
sented. An attempt was first made to understand the dynamics of a vortex/body interaction. It so happened that this 
is a rather complex process that led to some interesting properties, which were discussed in some detail. In this 
report only cases of zero angle of attack were investigated. However, extension to cases of nonzero angles of attack 
is relatively straightforward. 

The method of asymptotic matching enabled us to separate the problem into two parts: an inner potential flow 
region and an outer acoustics region. The main advantage of this is the efficacy of dealing with a complex geometry 
because only the potential flow sees the geometry, and there are efficient numerical methods for solving flow prob- 
lems of an arbitrary body. Matching was applied merely to transmit the gross properties of the inner region to the 
acoustic field. In fact the solution of acoustic pressure in the frequency domain is known; only the unknown coeffi- 
cients had to be determined by matching. 
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A fair amount of effort has been devoted to delineating the problems associated with the method of convolution 
integral and to demonstrate the usage of Hadamard’s definition of the final part. The evaluation of the convolution 
integral is not essential to solving acoustics problems, but it is of some interest in unraveling the vexing properties of 
a divergent integral. 

The evidence from examining calculated sound generation seems to indicate that sound radiation from a vortex/ 
bluff-body interaction is actually less than that from a corresponding vortex/airfoil interaction. Although this may 
appear to be paradoxical, there is an explanation for this property. We found a scaling factor that governs the sound 
generated by the vortex interaction and a mechanism to attenuate the background noise. In addition, we proposed a 
possible application that can lower sound production by simply positioning vortices on the “correct” side of the body 
instead of the “wrong” side. 
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Figure 1 .—Initial configurations of two, three, and four spinning vortices. 
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Figure 2. — Surface vorticity elements on a two- 
dimensional body. 
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Figure 3. — Comparison of surface velocity distributions of a potential flow over a step and 
over two rectangular bluff bodies. 
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Figure 4.— Trajectories of two spinning vortices passing an NACA 001 2 airfoil at zero angle of 
attack. For all cases, T = 0.1 5. Initial positions of vortices: (a) zq - -24.4 + /yo and z'q = -24.4 + 
/y‘ 0 ; (b) zq = -24.6 + /yo and z'q = -24.6 + iy' q; and (c) zo = -25.0 + /yo and z'o = -25.0 + iy l o, where 
yo = 0.29 and y’o = 0.1 2. 
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Figure 5. — Trajectories of two spinning vortices passing an NACA 0012 airfoil at zero angle of attack. 
Initial positions are the same as in figure 4(c): zq = -25.0 + 0.29/ and z'q = -25.0 + 0.1 2 /, except 
T = -0.15. 
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(a) 



Figure 6. — Trajectories of spinning vortices at zero angle of attack. For all cases, r = 0.15. (a) Two 
vortices passing an NACA 0012 airfoil. Initial positions: zq = -24.65 + 22/ and z'o = -24.65 + 0.12 /. 
(b) Three vortices passing a sharp-edged block. Initial positions: zq = -25.7 + 0.41/; z'o = -25.7 + 
0.24/, and z"q = -25.5528 + 0.325/. (c) Four spinning vortices passing an NACA 001 2 airfoil. Initial 
positions: zq = -25.72 + 0.41/; z'o = -25.72 + 0.24/; z"o = -25.55 + 0.24/; and z ,M o = -25.55 + 0.41/. 



x 


Figure 7. — Trajectories of two spinning vortices passing a sharp-edged block at zero angle of attack. 
Initial positions of vortices are the same for both cases: zq = -25.7 + 0.41/ and z'q = -25.7 + 0.24/. 

(a) T = 0.15. (b) T = -0.15. 
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Figure 8. — Trajectories of three spinning vortices at zero angle of attack. Initial conditions are identi- 
cal for both cases: T = 0.15; zq = -25.55 + 0.41/; z'q = -25.55 + 0.24/; and z" q = -25.4028 + 0.325 /. 
(a) Sharp-edged block, (b) Equivalent round-edged block. 



Figure 9. — Vortex trajectories rectified by a round-edged bluff body at zero angle of attack, initial 
conditions: f = 0.15; zq = -25.55 + 0.41/; z'o = -25.55 + 0.24/; and z" q = -25.19 + 0.325 /. 
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Figure 10. — Comparison of sound generation by two convected spinning vortices. {Mq = 0.2; dt = 0.025; 
point of observation: R = 9.0, 0 = 75° from positive x-axis.) (a) Analytic solution, (b) Method of con- 
volution integral. 
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Figure 1 1 . — Background noise due to a pair of spinning vortices attenuated by an NACA 001 2 airfoil 
(same as the background noise in fig. 13(c) but plotted on enlarged scale). 
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Figure 12. — Comparison of two methods for determining sound generation by a vortex inter- 
acting with Joukowski airfoil at zero angle of attack. Observer is directly above airfoil (6 = 90°). 
(a) Hardin and Mason’s results, (b) Results from discrete Fourier transform. 
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Figure 1 3. — Sound generation by a pair of spinning vortices interacting with an NACA 001 2 airfoil. (For all cases, 

T = 0.1 5, Mq = 0.2, dt = 0.01 25, and point of observation is R = 9.0, 0 = 75° from positive x-axis.) Initial positions 
of vortices are the same as in figures 4(a) to (c): (a) zg = -24.4 + iyo and z'g = -24.4 + iy‘ q; (b) zq = -24.6 + iyo and 
z'o = -24.6 + /y'o; and (c) z 0 = -25.0 + /y 0 and z' 0 = -25.0 + iy ‘ 0 , where y 0 = 0.29 and y' 0 = 0.1 2. 
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Figure 14. — Sound generation by a pair of spinning 
vortices interacting with an NACA 0012 airfoil. 

(For both cases Mq = 0.2; point of observation: 

R = 9.0, 0 = 75° from positive x-axis.) Initial positions 
are the same as in figure 4(c): zq = -25.0 + 0.29/ and 
z'o = -25.0 + 0.12/. (a) Convolution integral, r = 0.15, 
dt = 0.025. (b) Discrete Fourier transform, T = -0.15, 
dt = 0.01 25. 


NASATTM- 


1998-208403 


26 



0.6x10- 2 


0.4 b 


0.2 


p' 0.0 


- 0.2 


-0.4 


- 0.6 


— , 




J L 


. . w . 


5 10 16 20 25 30 


0.6x10-2 



Figure 1 5. — Sound generation by a pair of spinning 
vortices interacting with an NACA 0012 airfoil. 
(For both cases Mq = 0.2; dt = 0.0125; point of 
observation: R = 9.0.) (a) Faster spinning vortices 
with initial positions the same as in figure 6(a): 
zq = -24.65 + 22 / and z'q = -24.65 + 0.12/, 0 = 75° 
from positive x-axis. (b) Recalculation of sound 
generation for the same initial positions, except 
0 = 170°. 
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Figure 16. — Sound generation corresponding to tra- 
jectories in figure 7(a) for a pair of spinning vortices. 
(For both cases Mq = 0.2; point of observation: 

R = 9.0, 0 = 75° from positive x-axis.) (a) Discrete 
Fourier transform, dt = 0.01 25. (b) Convolution 
integral, dt = 0.025. 
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Figure 17. — Sound radiation by two spinning vortices 
interacting with rectangular blocks. Initial positions: 
T = 0.1 5, z 0 = -25.0 + 0.41/ and z* 0 = -25.0 + 0.24/. 
(For both cases Mq = 0.2; dt = 0.01 25; point of 
observation: R = 9.0, 0 = 75° from positive x-axis.) 
(a) Sharp-edged block, (b) Equivalent round-edged 
block. 
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Figure 1 8. — Sound generation by two spinning vorti- 
ces skirting a sharp-edged block. Initial positions: 

T = 0.1 5, zq = -25.0 + 0.52 / and z‘o = -25.0 + 0.35/. 
(Mq = 0.2; dt = 0.01 25; point of observation: R = 9.0, 
0 = 75° from positive x-axis.) 
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Figure 19. — Scaling factor for acoustic radiation from 
two and four spinning vortices interacting with an 
NACA 001 2 airfoil. (For both cases Mq = 0.2; 
off = 0.0125; point of observation: R = 9.0, 0 = 75° 
from positive x-axis.) (a) Two vortices and four 
vortices, (b) Four vortices corresponding to the 
initial positions in figure 6(c): T = 0.15, zq = -25.72 + 
0.41/; z'q = -25.72 + 0.24/; z" 0 = -25.55 + 0.24/; and 
z m 0 = -25.55 + 0.41/. 
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Figure 20. — Scaling factor for acoustic radiation from 
two and three spinning vortices interacting with a 
sharp-edged block. (For both cases Mq = 0.2; 
dt = 0.01 25; point of observation: R = 9.0, 0 = 75° 
from positive x-axis.) (a) Two vortices and three 
vortices, (b) Three vortices corresponding to the 
initial positions in figure 6(b): T = 0.15, zq = -25.7 + 
0.41/; z ' o = -25.7 + 0.24/, and z H 0 = -25.5528 + 0.325/. 
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